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The decay modes of the type B — > tttt are dynamically different. For the case 
B^ — > 7r + 7r~ there is a substantial factorized contribution which dominates. In 
contrast, the decay mode B® — > 7r°7r° has a small factorized contribution, being 
proportional to a small Wilson coefficient combination. However, for the decay 
mode B® — > there is a sizeable nonfactorizable (color suppressed) contribution 

due to soft (long distance) interactions, which dominate the amplitude. 

We estimate the branching ratio for the mode B^ — > 7r°7r° in the heavy quark 
limit for the b- quark. In order to estimate color suppressed contributions we treat 
the energetic light (u, d, s) quark within a variant of Large Energy Effective Theory 
combined with a recent extension of chiral quark models in terms of model- dependent 
gluon condensates. 

We find that our calculated color suppressed amplitude is suppressed by a factor of 
order h.QCD/mb with respect to the factorizable amplitude, as it should according to 
QCD-factorization. Further, for reasonable values of the constituent quark mass and 
the gluon condensate, the calculated nonfactorizable amplitude for B® — > 7r°7r° can 
easily accomodate the experimental value. Unfortunately, the color suppressed am- 
plitude is very sensitive to the values of these model dependent parameters. There- 
fore fine-tuning is necessary in order to obtain an amplitude compatible with the 
experimental result for B^ — > 7r 7r°. A possible link to the triangle anomaly is 
discussed. 
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I. INTRODUCTION 



Due to numerous experimental results coming from BaBar and Belle, there is presently 
great interest in decays of B-mesons. LHC will also provide us with more data for such 
processes. 5-decays of the type B — > titt and B — > Kir, where the energy release is big 
compared to the light meson masses (heavy to light transitions), has been treated within 
QCD factorization [1] and Soft Collinear Effective Theory (SCET) [2]. In the high energy 
limit, the amplitudes for such decay modes factorize into products of two matrix elements 
of weak currents, and some nonfactorizable corrections of order a s can be calculated per- 
turbatively. However, there are additional contributions of order AqcdA^£> which cannot 



be reliably calculated within perturbative theory [lL The so called pQCD-model and QCD 
sum rules have also been used for B- meson decays 3|, |4 ] . 

For decay modes which are of the heavy to heavy type, involving b- and c-quarks, the 
decay amplitudes have been described within Heavy Quark Effective Field Theory (HQEFT) 
0. Some transitions of heavy to heavy type in the heavy quark limits (l/m&) — > like B — B 
mixing |(| has been studied within Heavy Light Chiral Perturbation Theory (HL%PT) Q. 
Furthermore, other transitions which are formally heavy to heavy in the heavy quark limits 
(l/m&) — > and (I/77O — > 0, like the Isgur-Wise function 8| for B — >■ D, have been 
studied within HLxPT [7|. The cases B DD [9] and B D*^ [id, [ll| have also been 
studied within such a framework, even if the energy release in these processes is above the 
chiral symmetry breaking scale. Still this framework give amplitudes of the right order of 
magnitude. The calculation of such transitions have in addition been supplemented with 
calculations within a Heavy Light Chiral Quark Mod el ( HLyQM) to determine quantities 



which are not determined within HLxPT itse 
As pointed out in a series of papers 
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12|. 



13] , there are processes which have factorized 



amplitudes multiplied by a very small Wilson coefficient combination, such that nonfactor- 
ized amplitudes are expected to dominate. Examples are B® S ^-D°D° , B° — > D° rf 
12] and B® — > D°tt . The latter process B^ — > D° tt° was considered recently [3, 



In that case a heavy 6-quark decaying to a light, but energetic quark was involved. Then 



the light energetic quark might be described by an effective 
such a framework was Large Energy Effective Theory (LEET) 



;heqry. The first version of 



15 



16|. The HQEFT covers 



processes where the heavy quarks carry the main part of the momentum in each hadron. 
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To describe processes where energetic light quarks emerge from decays of heavy 6-quarks, 
LEET „ a8 Educed Q and U8 ed to 8tudy the current for r fl. 

The idea was that LEET should do for energetic light quarks what HQEFT did for heavy 
quarks. In HQEFT one splits off the heavy motion from the full heavy quark field, thus 
obtaining a reduced field depending on the velocity of the heavy quark. Similarly, in LEET 
one splits off the large energy from the full field of the energetic light quark, thus obtaining 
an effective description for a reduced light quark which depends on a light-like four vector. 
It was later shown that LEET in its initial formulation was incomplete and did not fully 
reproduce QCD physics Then LEET was further developed to be fully consistent with 
QCD and became the Soft Collinear Effective Theory (SCET) 

In the present paper we consider decay modes of the type B — > n it. The decay mode 
— > 7r~ 7r + has a substantial factorized amplitude, given by the current matrix element 
for B® — > 7r + transition times the matrix element of the weak current for the outgoing 
7r~, which is proportional to the pion decay constant The relevant Wilson coefficient 
is also the maximum possible, namely of order one times the relevant Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing factors and the Fermi coupling constant. This is in contrast 
to the process B® — > 7r°7r° which is color suppressed. As said above, decays of the type 
B — > 2ir have been extensively studied within QCD-factorization, SCET, and QCD sum 
rule methods [3]. In spite of tremendous efforts it has not been possible to obtain an 
amplitude compatible with the experimental result B® — > The purpose of this paper 

is study this decay mode within an alternative model dependent framework. 

First we point out that the factorized contribution to the decay mode B^ — > 7r°7r°, which 
is given by the B — > n transition amplitude times the decay constant of the 7r° meson, 
is almost zero because it is proportional to a very small Wilson coefficient combination. 
For the dominant nonfactorizable (color suppressed) amplitude for B^ — > 7r°7r° we will, as 
mentioned above, use a model named Large Energy Chiral Quark Model (LE^QM) recently 
constructed and used to handle the process B^ — > tt°D° Q- Here a variant of LEET 
was combined with ideas from previous chiral quark model (yQM) calculations similarly to 
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ill, Q 



19. 



what has been done for other nonleptonic decays 

A priori it might look strange to use the framework of chiral quark models when the energy 
release is big compared to the chiral symmetry breaking scale A x . The point is that the 
motion of the heavy quark or energetic light quark can be split off, and the various versions 
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of heavy-light or large energy chiral quark models and a corresponding chiral perturbation 
theory (x-PT) can be used to describe the redundant strong interactions corresponding to 
momenta of order 1 GeV and below. 

It might be argued that we should have used the full SCET theory as the basis our 
However, the purpose of our paper is to estimate, in analogy with previous 



new mode 
papers 



12 



19 



231 ]. the effects of soft gluon emission in terms of gluon condensates, 



where transverse quark momenta and collinear gluons will not play an essential role. In 
any case this construction [131] will be a model. Therefore it suffices for our purpose to use 



the more simp 



(XQM) |2lJ 



24 



e formulation of LEET. We will combine LEET with chiral quark models 



4271] . containing only soft gluons making condensates. In LExQM [13J an 



energetic quark is bound to a soft quark with an a priori unknown coupling, as proposed 



in 



211 ] . The unknown coupling is determined by calculating the known B — > n current 



matrix element within the model [13| . This fixes the unknown coupling because the matrix 
element of this current is known Then, in the next step, we use this coupling to 

calculate the nonfactorized (color suppressed) amplitude contribution to B% — > 7r°7r° in 



terms of t 
decay, fl 



lie 
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owest dimension gluon condensate, as have been done for other nonleptonic 
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19 



201 ] . After the quarks have been integrated out, we obtain an effective 



theory containing soft light mesons as in HL^PT, but also fields describing energetic light 
mesons. A similar idea with a combination of SCET with HL;\PT is considered in 



The LE^QM was constructed in analogy with the previous Heavy Light Chiral Quark Model 
(HLxQM) 20j and may be considered to be an extension of that model. 

One might think that to be completely consistent, we should also have calculated the 
Wilson coefficients within a relevant large energy framework. For this purpose the use of 
LEET would be dubious because it is an incomplete theory as mentioned above. However, 
as we will see below, the main uncertainty in our final amplitude will be due to uncertainty 
in our model dependent gluon condensate due to emission of soft gluons. Therefore the 
Wilson coefficients calculated within full QCD as in 29| will be appropriate for our purpose. 

In the next section (II) we present the weak four quark Lagrangian and its factorized 
and nonfactorizable matrix elements. In section III we present our version of LEET, and in 
section IV we present the new model LExQM to include energetic light quarks and mesons. 
In section V we calculate the nonfactorizable matrix elements due to soft gluons expressed 
through the (model dependent) quark condensate. In section VI we give the results and 
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conclusion. 



II. THE EFFECTIVE LAGRANGIAN AT QUARK LEVEL 

We will study decays of B% generated by the weak quark process b — > uud. We restrict 
ourselves to processes where the 6-quark decays. This means the quark level processes 
b — > duu . Processes where the anti- 6-quark decays proceed analogously. The effective weak 
Lagrangian at quark level is 29] (neglecting penguin operators) 



£ C ff = --^KbKTa I c a Qa + c B Qb] , 



(i) 



where the subscript L denotes the left-handed fields: q^ = Lq, where L = (1 — 75) /2 is the 
left-handed projector in Dirac-space. The local operator products Qa,b are defined as 



(2) 



In these operators summation over color is implied. In Eq. ([T]), ca and eg are Wilson 
coefficients. At tree level ca — 1 and cb = 0. At one loop level, a contribution to cb 
is also generated, and ca is slightly increased. These effects are handled in terms of the 
0tl0n Cr u P E^on S (RGE)fl, and the coefficient s can be calculated at for 
instance fi = or fi= 1 GeV. Using the color matrix identity 



2 ti n tfj — 5ij5in J-j-3in5lj , 

and Fierz rearrangement, the amplitudes for the processes B° — > 7r + 7r~ may be written as 



Gp 1 



4 —^V nh V* d 



ca + w c Cb ) ( n ~\ d Ll»u L \0)(n + \u L ^b L \B°) 
+2c B (n+n-\d Llfl t a u L u Ll n a b L \B )] , 



(3) 



and for B° — > 71 



°7l° 



G 



c B + jr c C A ) {AuLl^u L \Q){AdLYb L \B G ) 
+2 c A (vr 7r | d LlfM t a b L u Ll n a u L \B )] . 



(4) 



Here the terms proportional to 2ca and 2cb with color matrices inside the matrix elements 
are the genuinely nonfactorizable contributions. 



Since ca is of order one and cb of order —1/3 [12|, ll3j , we refer to the coefficients 

c f = (c A + -^c B ^J ~ 1.1 ; c n/ = + -j^ca^J ~ , (5) 

as favorable (c/) and nonfavorable (c„/) coefficients, respectively. Thus, the decay mode 
— * 7r+7r_ has a sizeable factorized amplitude proportional to c/. In contrast, the decay 
mode B^ — > 7r°7r° has a factorized amplitude proportional to the nonfavorable coefficient 
c n f which is close to zero. In this case we expect the nonfactorizable term (involving colour 
matrices) proportional to 2ca to be dominant, i.e. the last line of eq. (jl]) dominates. A sub- 
stantial part of this paper is dedicated to the calculation of this nonfactorizable contribution 
to the B® — » 7r°7r° decay amplitude. 

Thus the main task of this paper will be to calculate the matrix element of the operator 
Qc consisting of the product of two colored currents occurring in the last line of eq. (j4j): 

Qc = {d Ll ,t a b L ) {u Ll n a u L ) (6) 

for the color suppressed process B^ — > 7r°7r°. This matrix element will be estimated in 
section V where we use the LExQM to estimate nonfactorizable amplitudes in terms of 
emission of soft gluons making gluon condensates. 



III. AN ENERGETIC LIGHT QUARK EFFECTIVE DESCRIPTION (LEETtf) 

An energetic light quark might, similarly to a heavy quark, carry practically all the energy 
E of the meson it is a part of. The difference is that now the mass of the energetic quark is 
close to zero compared to the heavy quark mass rug and E, which are assumed to be of the 
same order of magnitude. We assume that the energetic light quark is emerging from the 
decay of a heavy quark Q with momentum pg = mq v + k. The heavy quark is described 
by the HQEFT Lagrangian for the reduced quark field Q v [5|]: 

£hqet = Qv {iv ■ D) Q v + 0{l/m Q ) , (7) 

where Q v is the reduced heavy quark field (often named h v in the literature), v its four 
velocity and tuq the mass of the heavy quark. 

The momentum of the light energetic light quark q can be written 

p^En^ + k^ , |A^| < \En>*\ m q < E , (8) 



7 



where E, which is of order mg, is the energy of the energetic light quark, m q is the light 
quark mass. Further, n is the light-like four vector wich might be chosen to have the space 
part along the z-axis, = (1; 0,0,1), in the frame of the heavy quark where v = (1,0). 
Then (v ■ n) = 1 and n 2 = 0. Inserting this in the regular quark propagator, In the limit 
where the approximations in (jBJ) are valid, we obtain the propagator 



S{p q 



7 ■ p q + m q 



7 ■ n 
2n-k 



(9) 



This propagator is the starting point for the Large Effective Theory (LEET) constructed in 



Ref. 



16] 



Unfortunately, the combination of LEET with xQM will lead to infrared di verg ent loop 



integrals for n 2 = (see section lIVp . Therefore, the formalism was modified 13|, |14| and 



instead of n 2 = 0, we use n 2 = 5 2 , with 5 = vjE where v ~ ^-qcd, such that 5 <C 1. An 
expansion in 5 will then within our model be equivalent to an expansion in Agc^/m^,. 
In the following we describe the modified LEET [16] where we keep 5^0 with (5 < 1. 



We call this construction LEET5 



13] and define the almost light -like vectors 



n = (1,0,0, +tj), ■ n= (1,0, 0,-77), 



(10) 



where rj = yl — 5 2 . This means that 



n 



M + = 2v tl , n 2 = n 2 = 5 2 , v ■ n = v ■ n = 1 , n ■ n = 2 - 5 2 



(11) 



In the following we use the projection operators given by 

V+ = -j^j ■ ■ n + 5) , V- = ■ n - 5)j ■ n 



(12) 



where N = \/2n ■ n = 2 + 0(5 2 ) . One factors out the main energy dependence, just as 
was analogously done in HQEFT, and define the projected reduced quark fields fl6| 



q±(x) = e^™V ± q{x) , q(x) = [q + (x) + q^x)} . 



(13) 



grangian 



13|: 



As in 16 1, the field g_ was eliminated and one obtained for q + = q n the effective La- 



Cleets <ln ( 7 j (in ■ D)q n + -^q n Xq n + 0(E 2 ) 



(14) 
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which (for 8 = 0) is the first part of the SCET Lagrangian. The operator X is given in 13 ]. 
Equation (Til]) yields the LEET5 quark propagator 



S n (k) = V, 



7 ■ n + 5 
N ' 



n ■ k) 



7 • n 
N(n ■ k) 



(15) 



which reduces to ([9]) in the limit 5 — > 0. In addition, for a light energetic quark, the 
propagator within SCET ^ will for small transverse quark momenta p± — > coincide with 
Eq. (JED. 

Based on LEET, it was found [16J in the formal limits Mjj — > oo and E — > oo, that 
a heavy H = (B, D) meson decaying by the weak hadronic vector current to a light 
pseudoscalar meson is described by a matrix element (P \ V 1 \ H) of the form 



{P\V\H) = 2E C (v) {M H ,E)n^ + ([ v> (M H , E) 



(«)/ 



where 



<"» = C^tf , C ~ (Aqcd) 3 ^ 2 



C (v) 



1 



(16) 



(17) 



This behaviour is consistent with the energetic quark having x close to one, where x is the 
quark momentum fraction of the outgoing pion 16]. 



IV. EXTENDED CHIRAL QUARK MODEL FOR HEAVY AND ENERGETIC 

LIGHT QUARKS (LExQM) 



The chiral quark model (xQM) 24], |25| and the Heavy-Light Chiral Quark Model 
(HL^QM) [20| , include meson-quark couplings and thereby allow us to calculate amplitudes 
and chiral Lagrangians for processes involving heavy quarks and low energy light quarks. In 
this section we will extend these models to include also hard, energetic light quarks. 

For the pure light and soft sector the xQM Lagrangian can be written as 19|, |24j: 



C 



xQM 



X [7 • (iD + V) + 7 • A - m) x 



(18) 



where m is the constituent mass term being due to chiral symmetry breaking. The small 
current mass term is neglected here. Here we have introduced the flavor rotated fields Xl,r- 



Xl = & qi 



XR = CQR , 



(19) 
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where q is the light quark flavor triplet and: 



V2 ^ n/6 



£ = exp{ffl//} 



n 



7T 



1T 



V 



7T 



7t_ j 2_ 



(20) 



Further, V M and ./L are vector and axial vector fields, given by 



(21) 



To couple the heavy quarks to mesons there are additional meson-quark couplings within 



HL X QM 



20|: 



Ant = -G H [Xa K Qv + Qv K Xa] , (22) 

where Q v is the (reduced) heavy quark field and H is the heavy (0", 1") meson field(s) 

#<+> =P+(v)( 1 -P* -H5P5) , (23) 

P* being the 1~ and P5 the 0~ fields, and P+(v) = (1 + 7 • w )/2. The quark-meson coupling 
Gh is determined within the HLxQM to be 20[ 



r 2 - 



;i + 3^) 



fy ' r 4(i+^f-^(^)) 

where 7?# = (8 — 7r)/64. The quantity p is of order one. 



(24) 



For hard light quarks and chiral quarks coupling to a hard light meson multiplet field 
M, one extends the ideas of xQM and HLxQM, and assume that the energetic light mesons 
couple to light quarks with a derivative coupling to an axial current |l3| : 



intg 



q W {id"M)q . 



(25) 



One combines LEET<5 with the xQM and assume that the ingoing light quark and the 
outgoing meson are energetic and have the behaviour exp (±iEn ■ x) as in ffl3|) . To describe 
(outgoing) light energetic mesons, we use an octet 3x3 matrix field M = exp (+iEn ■ x) M n 
of the same form as II in (1201): 



(26) 





( 1 Vn 


7T + 

n 


K+ \ 

n 


M n = 




"J 1 Vn 


K° 

n 




K K- 


n 


V6 J 
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where 7r°, 7r+, K£ etc. are the energetic light meson fields with momentum ~ En 11 . 

Combining f l25|) with the use of the rotated soft quark fields in ffl9|) and using <9 M — > iE n M 
one arrives at the ansatz for the LE^QM interaction Lagrangian: 

C-intqS = G A E\ {j-n)Zq n + h.c. , (27) 

where q n represents an energetic light quark having momentum fraction close to one and x 
represents a soft quark (see Eq. ffl9|) ). Further, the coupling Ga is determined by physical 



requirements 



13 



16j . and 

Z = £M R R-£}M L L . (28) 



Here and Mr are both equal to M n , but they have formally different transformation 
properties, This is analogous to the use of quark mass matrices Ai q and in standard 
Chiral Perturbation Theory (%PT). They are in practice equal, but have formally different 
transformation properties. 

The axial vector coupling introduces a factor 7-71 to the vertex (see (|27|) ). which simplifies 
the Dirac algebra within the loop integrals. In order to calculate the nonfactorizable contri- 
bution, one must first find a value for the large energy light quark bozonisation coupling Ga- 
This was done [l3] by requiring that our model should be consistent with the equations ( fl6|) 
and (EETj). Applying the Feynman rules of LExQM [l3| we obtain the following bosonized 
current (before soft gluon emission forming a condensate is taken into account): 

Jo(H Vb -> M n ) = -N C J dkTr{YLiS v (k) [-iG H H^] iS x {k) [iE Ga 7 -nZ) iS n (k) }(29) 

where dk = d D k/(2%) D (D being the dimension of space-time), and 



are the propagators for heavy quarks described by (jlBj) . for light constituent quarks, and 
(I14p for light energetic quarks. The presence of the left projection operator L in Z ensures 
that we only get contributions from the left-handed part of the interaction in f )27|) . that is, 
Z — > —£*Ml L. The contribution in ( )29|) corresponding to the B — > 71 current is illustrated 
by the lower part of the diagram in Fig. [TJ 

Loop diagrams within LExQM depend on momentum integrals of the form 

r dk 
rst = J (v ■ ky (k ■ n y (f - m?y ' ( } 
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FIG. 1: The factorized contribution to the B — > tt + it~ decay, as described in combined HLxQM 
and LExQM. Double lines, single lines and the single line with two arrows are representing heavy 
quarks, light soft quarks and light energetic quarks, respectively. Heavy mesons are represented 
by a sigle line combined with a parallel dashed line, and a light energetic pion is represented by a 
dashed line with double arrow. 



rst ~J (vky(k-nY(k 2 -m 2 y~ rst + rst ' 1 ] 

These integrals have the important property that K^} dominates over and K rst with 
one power of 1/5. In the present model, we choose v = m which is of order Aqcd- Thus the 
constituent light quark mass m is the equaivalent of Aq C d within our model. Some details 



13|. 



of the calculation of the B — > n is given in Ref. 

To calculate emission of soft gluons we have used the framework of Novikov et al. [30]. 
In that framework the ordinary vertex containing the gluon field A ^ will be replaced by the 
soft-gluon version containing the soft gluon field tensor G a 



igjrVAl^ -\gst a T^G a ^....\k=o , (33) 



where k is the momentum of the soft gluon. (Using this framework one has to be careful 
with the momentum routing because the gauge where A ^ = has been used.) Here 
T^ = 7^ , v 1 *, or (7 ■ h + S)/N for a light soft quark, heavy quark, or light energetic 
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quark, respectively. Our loop integrals are a priori depending on the gluon momenta ki 2 
which are sitting in some propagators. These gluon momenta disappear after having used 
the procedure in ( 1331) . (Note that the derivative has to be taken with respect to the whole 
loop integral). 

Emision from the heavy quark or light energetic quark are expected to be suppressed. 
This will be realized in most cases because the gluon tensor is antisymmetric, and therefore 
such contributions are often proportional to 

G a ta ,v ti v v = , or GJXn" = 0. (34) 

However, there are also contributions proportional to : 

G%,v*n v ? . (35) 

analogous to what happens in some diagrams for the Isgur-Wise diagram where there are 



two different velocities vj, and v c 3JJ. Such contributions appear within our calculation when 



two soft gluons are emitted from the heavy quark line. 



Using the prescription 19 



2Q|, 



25 



30) 



9 2 s G%G a pX -> 4vr 2 (^G 2 )^(^ A - g^ 9l/p ), (36) 

for the gluon condensate one obtains the leading bosonized current [1 

JUH -> M) = -i^^m 2 FTr {^LH^ [ 7 • n]?M L ] , (37) 

where the quantity F obtained from loop integration is a priori containing a linearly divergent 
integral, which is related to the axial coupling g^, and can be traded for g^. One obtains 



131 ] for the quantity F: 



3f« n v N c (24-7tt) .a a _2. 



, > v 1 "^) + T7T " 7fi « 4 - (-° ) ■ (38) 
8m^p 167T 768 m 4 7r 

Note that F is dimensionless. The parameter p is given in ff24l . Numerically, it was found 
3 that F ~ 0.08. 

In order to obtain the HLxPT Lagrangian terms Tr(H a H b v^V^ a ) and Tr(H a H b/ ~f f j i , ~f 5 A^ a ), 
having coefficients +1 and —g^ respectively, one calculates quark loops with attached heavy 
meson fields and vector and axial vector fields V 1 or A^ 1 . Then logarithmic and linearly 



13 



divergent integrals obtained within the loop diagrams are identified with physical quantities 



or quantities of the model 



19 



20 



24 



25|. 



In order to fix Ga in ( 12T1) . we compare ( fl6l) with ( 1371) . In our case where no extra soft 
pions are going out, we put £ — > 1, and for the momentum space Ml — > ku^B-, with 
the isospin factor ku = 1/ v^2 for 7r° (while ku = 1 for charged pions) . Moreover for the 
-B-meson with spin-parity 0~ we have Hv — > P+ (v)(— i'ysj^Mff. Using this, the involved 
traces are easily calculated, and we obtain Jt ot {H — > M) for the case B® — > 7r + : 

JL<M^* + ) = ^^{^/M h 7E)m'Fn^. (39) 



Using the equations (fT6l) . fl38l) . and (1391) . one obtains [13] 



n 

where £^ is numerically known [32]. Within our model, the analogue of Aqcd is the 
constituent light quark mass m. To see the behaviour of Ga in terms of the energy E, the 
quantity C in (TT71) is written as C = cm 2, which gives 

GU = ( %= 1 4 , (41) 



which explicitly displays the behaviour ~ E~ 3 / 2 . In terms of the number iV c of colors, 
A ~ v^Vc and F ~ N c which gives the behaviour Ga ~ 1/ y/Nc, i.e. the same behaviour as 
the coupling Gh in (1221) . 

The bosonized current in (1371) can now be written as 



JUH M) = -2<^ y/ A Tr { 7 ^ +) [ 7 • n] ?M L } . (42) 

V. NONFACTORIZABLE PROCESSE IN LExQM 

In this section we calculate the nonfactorizable contribution to B® — > 7r°7r° in Eq. (jlj). 
This will be formulated as a quasi-factorized product of two coloured currents, as illustrated 
in Fig. [2j Then the nonfactorized aspects enters through color correlation between the two 
oarts, using Eq. ( 1361) . Such a calculation within HL^QM and HL%PT is done previously 



9j for B® s — > D°D°. Here we will use the colored current for B — > tt, within the LExQM 
presented in the preceding section; see the diagram in Fig. [21 Using the Ga value from the 
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7T 



FIG. 2: Nonfactorizable contribution containing large energy light fermions and mesons. There is 
also corresponding diagram where the outgoing anti-quark u is hard. 



preceeding section, we may now calculate the nonfactorizable contribution to the process by 
adding one soft gluon to each loop. Then we calculate the decay width for B® — > 7r°7r° from 
this nonfactorizable amplitude, and compare our results with experiment. 



For a low energy guar 
effective propagator 
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i interacting with one soft gluon, one might in simple cases use the 



33] 



S?(k) = ^-t a G 



g. aria (2ma^ +{a^, 1 -k}) 



(k 2 



m 



2\2 



(43) 



where {a, b} = ab + ba denotes the anticommutator. This expression is constistent with the 
, and can be used for the diagram in Fig. |2] . 



prescription in (133 

Then one gets 13| the following contribution to the bosonized colored B — > tt current, 



shown in the lower part of the diagram in Fig. |2] : 



J? G (H -> M) a 



dkTr{YLt a iS v {k) [-iG H H^] iS?(k) [%EG Al -nZ\ iS n {k) }(44) 



where a is a color octet index. Once more, we deal with the momentum integrals of the 
types in ( 13T1) and (|32l) . Taking the color trace, rewriting OHj), we obtain a contribution of 
the form 



J? G (H b M) a = g s G a a0 T^(H b M) 



(45) 
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where the contribution from the (lower part of) the diagram in Fig. [2] : alone is to leading 
order in 5 

T^(H b -)■ M) = %^^ A n CT Tr {-fLH™ lx t*M L ) , (46) 

where E ■ 5 = m has been explicitly used. 

There are also other diagrams not shown. In one case the gluon is emitted from the 
energetic quark. This diagram is zero due to (1341) . Furthermore, there is a diagram not 
shown where the gluon is emitted from the heavy quark which contains a non-zero part 
due to (1351) . This gives an additional contribution to the colored B — > n current which is 
nonzero. However, this one will be projected out because it should be proportional to the 
Levi-Civita tensor to give a nonzero result for the B® — > 7r°7r° amplitude as a whole, as will 
be seen from Eq. (j4"9~|) below. 

The colored current for an outgoing tt° should now be calculated in the LExQM (see 
upper part of the diagram in Fig. [2]), and we find 

Jf G (M fl )° = - J dkTr {YLt a iSf(k) [i EG a 7 • n Z\ iS h (k) } , (47) 

This colored ir° current has the general form 

J?c(M h ) a = g s G a af3 T^(M h ) , (48) 

where the tensor T is given by 

T^(M n ) = 2 Yn a e aa ^ Tr [\ x M h ] , (49) 

where the X x within the trace is the appropriate Gell-Mann SU(3) flavor matrix. For an 
outgoing hard tt° this trace has the value \JE~]2 when going to the momentum space. The 
explicit factor 2 in front of this expression comes from the corresponding diagram, where 
in the upper part of the diagram the antiquark could be hard and the quark could be soft 
and emit a soft gluon. The factor Y contains the result of loop momentum integration. The 
relevant loop integral is now 



which gives 



_ dkk» _i 2 
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Here the parameter A is of order 10~ 2 to 10 _1 and very sensitive to small variations in the 
model dependent parameters m and (^-G 2 ). 

It is easily seen that the experimental value of the B® — > 7r°7r° amplitude can be accomo- 
dated for a constituent mass m around 220 Mev and a value for (^G 2 )i /4 around 315 MeV. 
These values are of the same order as used in previous articles p, O, 11-13. 20-22]. But in 



contrast to these previous cases the present amplitude for B® — > 7r°7r° is very sensitive to 
variations of the model dependent parameters m and (^-G 2 ). Or more specific, the colored 
current J^ G (Mn) a in ( |4"T|) . (HHj) . and ( H9|) is very sensitive to these parameters. In other 
words, Y\ has to be fine-tuned in order to produce the experimental result. 

In a recent paper [3J| an extra mass parameter was introduced in the propagator of heavy 
quarks. One might do the same for propagator of the light energetic quark, and use 

Sn = N(n \ + A B ) ' (52) 

This would also bring this propagator more in harmony with the SCET propagator if A n ~ 
p\/E . This will to first order in A n give an extra contribution in the loop integral obtained 
from the diagram in Fig. 2: However, taking into account also the corresponding diagram 
where the light anti-quark is the energetic one, this first order term in A„ cancels. But there 
will be terms of second order in A n , which are of order 5 2 . Such contributions have to be 
considered together with higher order (in 5) terms obtained from the interaction given by 
the operator X in (j!4p . 

One should note that the colored current given by (148]) and (1491) is determined by a 
triangle diagram. Thus one might speculate if it can in some way be related to the triangle 
anomaly. Namely, the diagram in Fig. [2] would have, for standard full propagators, the 
mathematical properties of the diagram relevant for the triangle anomaly. Using dimensional 
regularization in this case, with dimension D = 4 — 2e, the loop integration gives an divergent 
result ~ I2 ~ 1/e while the corresponding Dirac trace is ~ e. Thereby one obtains a finite 
expression for the triangle diagram in that case. However, in the present case we have 
replaced one of the standard (full) quark propagators by the SCET-like propagator S^. 
Then the trace will not be ~ e while the corresponding loop integral is still divergent ~ 1/e. 
This means that the diagram is in total divergent. Within our various chiral quark models 
including heavy quarks and light energetic quarks, the naive dimensional regularization 
(NDR) has been used, and divergent integrals have been identified with physical parameters 
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11 



13 



20 



22 



31] . Using other schemes additional finite terms of type e/e might appear 



191 ]. and some parameters might have to be redefined. 
We also note that the description of the anomaly is rather tricky when going from the 



ow energy process tt° — > 27 to higher energies where some cancellations occurr 



35j. In 



34j the high energy processes Z — > 7r°7 and 7* — > 7r°7 was studied. (Here the high energy 



virtual photon 7* is coming from an energetic e + e~ pair). In this case a part of the amplitude 
corresponding to low-energy decay 7r° — > 27 is cancelled. But there is a remaining anomaly 
tail relevant for some high energy processes 34j, |35j . Trying to adapt such a description in 



our case, the tensor T in ( 1491) for an outgoing 7r° and soft gluon would be replaced by 

lAn 



T" ap {An) 



7T aa/3a 
fa c J 



(53) 



where we have taken into account that couplings and color traces are different from the 



calculations in 



34 



35]. The quantity I An is an integral given by 

oc doc 



~ > (54) 

7]X{1 — X) — 1 

where 77 = p 2 /m 2 . Using as before m as a constituent mass and p 71 " = Eh would give ri = 1 



leading to I An — 0.6. However, as the anomaly tail is of perturbative character |34j . |35J one 
might think that it is more relevant to use masses closer to the current masses of order 5 to 
10 MeV. In this case one has an assymtotic behaviour J^ n ~ ln(r])/r] , and this would give 
values for Ia u of order 10~ 2 . 

Now we use ( 1361) and also include the Fermi coupling the Cabibbo-Kobayashi-Maskawa 
matrix elements, and the coefficient 2ca for the nonfactorizable contributions to the ampli- 
tude, where C4 is the Wilson coefficient for the Oa local operator. Using Eqs. (T45|) and 
ff4T|) we find the effective Lagrangian at mesonic level for the nonfactorizable contribution 



to B° d 7r°7r° 



r LE x QM 
^Non.fact. 



47T 2 C4 



^v ub v: d 



'—G 2 ) S(H b — >■ M n Mn) 

71 



where S(H b — >■ M n M n ) is the tensor product 

S(H b -+M n M h ) = T^(H b ^M n )T^(M n 



(55) 



(56) 



Using Eqs. (T46]) and P9l) . and n-h ~ 2 , we find the amplitude expressed entirely by known 
parameters, we find an explicite expression for S(H b — > M Ma) in the case — > 7r°7r : 

1 2 G\ Gh 



6 



y/2 1287T 



Y E 



M. 



B ■ 



(57) 
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We will now compare this nonfactorizable amplitude for B®, — > 7r°7r° with the factorized 



amplitude which dominates B® — > 7r + 7r 



AW- = ( 4 ^^* d ) • c f ■ Q j M (o) • Q^°^ + )) 



(58) 



where 

J^tt") = /.M, , J"CB°->7r + ) = 2£r^C (,;) • (59) 



The form factor Q v > is defined in f fl6|) and (JT 

Using the equations (HUj) and ( 155]) - ( 159|) . we find the following ratio between the non 
factorized for B® — > 7r 7r and the factorized amplitudes B^ — > 7r + 7r~ is 



M(£?° 7r+7r-) Fact c/ iV c v^M^ 

where k is a model-dependent hadronic factor 



(60) 



ttAU^G 2 )Y , 
re = V", 7= • (61) 



It will be interesting how the ratio r scales with energy E. Using the scaling behaviour for 
with C = cmh in (fTTI) we find for the ratio r: 

Ca kc m , . 

Our calculations show that the ratio r of the amplitudes are suppresed by 1/N C , as it should. 
The ratio is also scaling like m/E. Because E ~ rrib/2 and m is the equivalent of Aqcd in 
our model, we have found that the nonfactorized amplitude is suppressed by Aq C D/ m b as 
required by the analysis in ref. jj. 

Concerning numerical predictions from our model, we have to stick to Eq. fl60l . The 
measured branching ratios for B2—> 7r~7r + and B® — > 7r°7r° are (5.13 ± 0.24) x 10 -6 and 
(1.62 ±0.31) x 10 -6 , respectively |36| . In order to predict the experimental value solely with 
the mechanism considered in this section, we should have r ~ 0.56 ± 0.11. Numerically, we 
use ~ 1/3 [32j |. In previous papers on the heavy-light chiral quark model constituent 
masses m ~ 220 MeV and (^G 2 ) 1 ^ 4 ~ 315 MeV has been used. From the plot of r in 
Fig. [3j we observe that the experimental value of r can easily be accomodated by values of 
such orders. The bad news is that in our case the value of Y\ and thereby k and r is very 
sensitive to the explicit choice of m and (^-G 2 ) 1 ^. Thus fine-tuning has to be used. 
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FIG. 3: Plots for the quantities F, p, A and r in terms of m and (^G 2 ) 1 . We observe that 
for reasonable values of these parameters the ratio r can take a wide range of values such that 
fine-tuning is required to reproduce the experimental value. 



We also find that the perturbative anomaly tail will numerically reproduce the amplitude 
for I An — 3.2 x 10~ 2 , corresponding to a quark mass mo — 11 MeV, i.e. of same order of 
magnitude as typical current quark masses. Using a hybrid description with a quark model 
with constituent quark masses for the colored B° — > 7r° current in (H4l)-( H6l) . and the anomaly 



tail description 34j, |35] for the colored n° current in (H71) . ()48l and ( )49l) is not preferrable. 



Also, such a hybrid description also fails to show the behaviour Ag C D/ m b required by QCD- 
factorization. Still it might be interesting that we can numerically match the colored current 
for outgoing tt° with the anomaly tail description. . 

Note that there are also mesonic loop contributions similar to those contributing to 



processes of the type B — )■ D D and B — > 7 D [9_|, 111] . For those processes intermediate 
D*(l~) mesons contributed. In the present case the analogous contributions would involve 
energetic vector mesons p n , and we would need the amplitudes for B — > p n p n . Such loops 
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FIG. 4: Meson loops for — > tttt. The zig-zag lines represent energetic p-mesons. The dashed 
lines with double arrow are energetic light mesons and the dashed line with no arrow is a soft pion. 



are shown in Fig. |H The diagram to the right would be calculable within an extended theory 
involving energetic vector mesons. Unfortunately while the diagram to the left would be 
dubious because typical loop momenta would significantly exceed 1 GeV, and would require 
insertion of ad hoc form factors or should be handled within dispersion relation techniques. 
Both diagrams would of course require knowledge of the p n ir n n coupling in Fig. HI In any 
case such calculations are beyond the scope of this paper. 



VI. CONCLUSION 



We have pointed out that the factorized amplitude for process B® — > tt°it is proportional 
to a Wilson coefficient combination close to zero. Thus the nonfactorizable contributions 
dominate the amplitude for this decay mode. To handle the nonfactorizable contributions 



we have extended previous chiral quar 



i models for t 



m 
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25], and the heavy light case 



le pure 



20j used in 



q, y. 



ii 



light quark case[2J| used in 



12 
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221 ]. to include also 



energetic light quarks. 

We have found that within our model we can account for the amplitude needed to explain 
the experimental branching ratio for B® — > 7r°7r° [36|. In addition the ratio r between 
the non-factorizable and factorized amplitude scales as KQcoj^b in agreement with QCD 
factorization [1]. However, the bad news is that the calculated amplitude is very sensitive 
to our model-dependent parameters, i.e. the constituent quark mass m, and the gluon 
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condensate (^-G 2 ). Anyway, final state interactions should be present 37l |. 
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